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Hadronic Correlators from All-point Quark Propagators 
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A method for computing all-point quark propagators is applied to a variety of processes of physical interest in 
lattice QCD. The method allows, for example, efficient calculation of disconnected parts and full momentum-space 
2 and 3 point functions. Examples discussed include: extraction of chiral Lagrangian parameters from current 
correlators, the pion form factor, and the unquenched eta-prime. 



1. Methodology 

The very high cost of generating decorrelated 
dynamical gauge configurations makes it increas- 
ingly important to extract the maximum physi- 
cal information content of each available config- 
uration. In many cases, this requires calculation 
of hadronic observables involving quark propaga- 
tors from any point on the lattice to any other. 
Here we describe an approach to obtaining such 
propagators by simulating bosonic pseudofermion 
fields. Introduce bosonic pseudofermion field 4>„ia 
with action (m a lattice site, a the spin-color in- 
dex, Q the Wilson or clover operator): 
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For fixed background gauge field A, simulating 
the pseudofermion field with the preceding ac- 
tion produces the following correlator (<< O » 
means the average of O relative to the measure 
e-^) : 

<< 4>rna(pnb »S{4,) = {H^'^)ma,nb 
« <t)ma{4>^ H)nb » S(<j}) = {H ^)ma,nb 

^b)ma^nb 

These simulations are practical for two reasons: 
(i) The pseudofermion average << .... >> is 
efficiently implemented by heat-bath update of 
pseudofermion fields. 



(ii) For a fixed gauge field, most quantities decor- 
relate after a few pseudofermion sweeps. 
The computation of multipoint hadronic corre- 
lators involving n quark propagators can be re- 
duced to convolutions of n pseudofermion fields, 
rapidly computed by fast Fourier transform 
(FFT). For example, the full 4-momentum trans- 
form A((7) = Ex,ye*«'^''"^^A(x,7;) of the 2-point 
pseudoscalar correlator 

A(x,j/) = <0|r{*(a:)75^'(a:)*(y)75*(2/)}|0> 

= - < tr((Q-l75)xy(Q"S5)yx) > 

= - <Pxai<p''H)ybXybix''H)^a » 

ab 

= {(!)'' Hx)yy{x^H(t))^^ » 

becomes 

A{q) = - « FFT(x^H(/))(q)FFT(?itHx)(-q) » 

2. Applications 

We have studied the feasibility of a pseud- 
ofermion approach to all-point propagators in a 
number of examples of physical interest. A brief 
report on the progress to date follows. 

2.1. Chiral Lagrangian Parameters from 
Current Correlators 

The chiral Lagrangian predicts the low- 
momentum structure of 2-point functions such as 
A{q) in QCD 0: 
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configuration 10000 
configuration 46950 
average A{q^) 
fit {q^:1 -10, lattice units) 




Figure 1 . Decorrelation of pseudofermion evalua- 
tions of A{q) 



Figure 2. 4 parameter fit of measured A{q'^) 
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We have studied various correlators of this type 
using 800 dynamical configurations generated us- 
ing the truncated determinant (TDA) algorithm 
1^ on large coarse (6^) lattices. The presence of 
low eigenmodes of H leads to longer autocorre- 
lation times for the zero momentum component 
(see Fig.l), on the order of 60 pseudofermion 
sweeps, with much more rapid decorrelation for 
nonzero momentum values of A{q^). (Low mo- 
mentum autocorrelations can be substantially re- 
duced by projecting out the lowest few eigen- 
modes of H- the required code has been written 
and is presently being tested [^). However, even 
at zero momentum, the intrinsic gauge fluctua- 
tions exceed the statistical errors from the pseud- 
ofermion evaluation (which become very small at 
higher momenta where decorrelation is rapid). 
This can be seen in Fig. 2, where we display the 
measured A(g^) for two separate configurations, 
as well as the average over 800 configurations and 
a fit to the chiral formula. In fact a good fit to the 
predicted chiral form can be obtained even with- 
out the zero- momentum point. The parameters 
hijli are couplings in the next-to- leading order 
chiral Lagrangian Fitting the measured A(g) 
(800 6'' unquenched TDA lattices |) with 0(0^) 



improved gauge action, a=G.4 F): 

+ A2 

A standard cosh fit of smcarcd-local correlators 
gives A-U = 0.396 ± 0.007. Allowing the pion 
mass to vary in the chiral formula, the best fit 
(see Fig. 2) is obtained in the range 0.25< q^ <2.5 
GeV^, and gives a 2% evaluation (jackknife er- 
rors) of the one-loop chiral parameter ^3: 

Ai = 23.1 ±0.4 

A2 = 0.178 ± 0.017 =^ = 0.422 ±0.020 

^3 = 9.48 ±0.16 

Ai ^ -0.70 ±0.02, ^ = 0.025 ±0.001 

2.2. Three-point Functions: the Pion Form 
Factor 

To extract the pion formfactor, we need the 
following 3-point function (see Fig. 3): 

wxyz 

< §(z + f, t2)75*(x, i2)^(y, ii)7o^'(y, h) 
§(u;,io)75*(0,to) > 

To evaluate J{q^) simultaneously for all space- 
like q injected at time ti by the electromag- 
netic current, it suffices to use a conventional 
smeared source propagator for the quark prop- 
agation from w, to to y, ti and a smeared sink 
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Figure 3. 3-poiiit function for pion form factor 



propagator for the propagation from x + z,t2 to 
0, to- The remaining propagator from y,ti to x,t2 
is then needed for all source and sink points (to 
obtain the full momentum space Fourier trans- 
form) and is evaluated by the pseudofermion tech- 
nique. At this stage, both the connected and dis- 
connected contributions (where the quark prop- 
agates from ?7 back to y) to J{q^) are readily 
available, as the all-point propagator also gives us 
the amplitude for all coincident source and sink 
points. An approximate pion form factor is given 
by J(<72)e(*^-*i)^(9)/i;(g)) (to, ^i, ^2=0,3,6, with 
E{q) the energy for a lattice pion of momentum 
q): results obtained with 60 quenched 12''x24 lat- 
tices at (3 =5.9 are shown in Fig. 4. Analysis 
of a single gauge configuration takes about 12 
hours on a Pentium 4 processor. Reliable cal- 
culation of the pion form factor requires use of 
optimized smearing wavefunctions f^"^ to project 
out ground-state pions as the signal dies quickly 
at larger Euclidean times ti — to,t2 — h, espe- 
cially at larger momentum, and to check that 
a plateau is reached at large Euclidean time. 
Of course, more accurate results require higher 
statistics than used in this feasibility study. 

2.3. The unquenched eta-prime 

We may extract the ctaprimc propagator in the 
isoscalar channel using two pseudofermion fields 
to generate propagators for both quark lines: 
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Figure 4. Connected and Disconnected contribu- 
tions to Pion Form Factor 



In terms of pseudofermion averages, this gives a 

connected 

xyzw 

« (l>{x){^^ H){y + w)x{y){x} H){x + z) » 

as well as disconnected ("hairpin") contribution 
(for two flavors of sea quarks): 

-2 Y e'P<''-y^f'^{z)r"'{w) 
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« (i>{x){4>^ H){x + z)x{y){rH){y + w) » 

Computations of the unquenched eta-prime prop- 
agator using this method are in progress, using 
10^x20 lattices generated with the TDA algo- 
rithm. 
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